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Abstract 

In the history of condensed matter physics, reinvestigation of a well-studied material with en¬ 
hanced quality sometimes led to important scientific discoveries. A well-known example is the 
discovery of fractional quantum Hall effect in high quality GaAs/AlGaAs heterojunctions. Here 
we report the first single crystal growth and magnetoresistance (MR) measurements of the silver 
chalcogenide /3-Ag2Se (Naumannite), a compound has been known for the unusual, linear-field- 
dependent MR in its polycrystalline form for over a decade 1, 2]. With the quantum limit (QL) 
as low as 3 Tesla, a moderate held produced by a superconductor magnet available in many labo¬ 
ratories can easily drive the electrons in Ag 2 Se to an unprecedented state. We observed significant 
negative longitudinal MR beyond the QL, which was understood as a ‘charge-pumping’ effect be¬ 
tween the novel fermions with opposite chiralities. Characterization of the single-crystalline Ag 2 Se 
and the fabrication of electric devices working above the QL, will represent a new direction for the 
study of these exotic electrons. 
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Electrons in metal move in the quantized orbits known as Landau levels (LLs) in an exter¬ 
nal magnetic field [3, 4], The so-called quantumn limit(QL) is attained when the magnetic 
field is higher than that corresponds to the lowest LL, in which all the electrons are confined 
in one highly degenerate state. Beyond the QL, the electrons are further fractionalized in 
two-dimentional(2D) electron gas, while little is known about the three-dimentional(3D) sys¬ 
tems. Unfortunately very few materals manifest their electron concertration dilute enough 
so that the QL is accessble. Here we show that the single-crystalline /3-Ag2Se, a well-known 
quantum material, is one candidate for the exporation of this unknown region. 

The mineral Naumannite (/3-Ag 2 Se) was first discovered in 1828. It crystallizes in a 
non-central-symmetric structure which has an orthorhombic unit cell (P2 1 2 1 2 1 ) with two 
crystallographic distinct silver atoms and one selenium atom j3, H] (Fig. H] insets). Pre¬ 
vious studies showed that polycrystalline /3-Ag 2 Se manifested large, linear-field-dependent, 
transversal magnetoresistance (MR) at both low and high temperature regions 0,0. This 
linear MR remains unsaturated as a function of magnetic field strength up to 50 Tesla (T). It 
was believed that the origin of this unusual linear MR is either a quantum effect of electron 
correlation due to a linear dispersion at the energy where the conduction and valence bands 
cross (2, 8], or a classical effect of spatial conductivity fluctuations in strongly inhomoge¬ 
neous polycrystals (9, 10]. Very recently topological characteristics of electrons in Ag 2 Se due 
to spin orbital coupling (SOC) were emphasized jll |. 

Although various papers have been published on the studies of polycrystalline and nano¬ 
size Ag 2 Se 12-|l4j, research on a macro-size single crystal has never been reported until now. 
The difficulty of the single crystal growth comes from the thermal instability of the /3-phase: 
it undergoes a first order, polymorphic structural transition to a cubic, high-temperature a- 
pliase at 133° C. In the a-phase, selenium atoms form a body-centered cubic frame and silver 
atoms are distributed over several interstitial sites jl5, 16]. During a regular crystal growth 
process, ct-Ag 2 Se intends to form a multi-domain polycrystalline /3-phase when cooling to the 
room temperature. By contrast, we discovered a novel modified self-selecting vapour growth 
technique [i0, which enables us to obtain single crystals of /3-Ag 2 Se above the transition 
temperature (See details in Method Part and Supplemental Information (SI)). During this 
growth process, the vapour from polycrystalline Ag 2 Se condensed on small seed plates of 
the Q-phase, leading to the ribbon-like /3-phase single crystals (Fig. [I] inset). Macro-size 
single crystals of /3-Ag 2 Se exclusively provide a chance for studying the electronic states 
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via transport measurements with less influence of scattering from defects, particularly from 
the grain boundaries in polycrystals. We focus on the MR measurements for the current 
(/) along the crystallographic a direction while the magnetic field can be along different 
directions in this paper (Fig. |T] inset). 


Temperature dependent resistivity of sing 


that of previously reported in polycrystals 


e crystalline Ag 2 Se shows a similar profile as 
, 7]: a broad maximum occurs around 100 K, 


below which the resistivity shows a metallic-like behavior (Fig. [[]). The resistivity of the 
samples R1 and R2 that are from one growth batch shows one order of magnitude less than 
that of the samples R3 and R4 from the other. Hall measurements (for ///a, H// b, and 
Hall voltage along c) reveal that all samples are n-type for the whole temperature range 
and the carrier density is n — 4.9 x 10 19 cm -3 and 2.4 x 10 18 cm -3 for the samples R1 and 
R4 at 2 K, respectively (See SI and Table []])• The n-type carriers most likely come from a 
small amount of deficiencies of selenium atoms commonly observed in selenide compounds 
[lsj ]. Both R1 and R4 show mobility (fib = 1/pne) about 1000 cm 2 /(V • s) at 2 K, leading 
to a product of the classical cyclotron frequency to the scattering time, cocT = HbB < 1 at 
9 T [3|. 

Figure [2] shows the MR (defined as A Ph/Po) °f the samples R2 and R4 along the three 
principal axes at representative temperatures. The transversal MR for H //b and c shows 
a similar profile: it crosses over from a quadratic dependence on H at low fields, to a linear 
and unsaturated dependence up to 9 T. This linear and unsaturated MR up to 300 K is 


similar to the previous reports in polycrystals 


fl H, 


19]. For H //b in R1 and R2, weak 


Shubnikov-de Haas (SdH) oscillations occur on the linear background when T < 10 K and 
fi 0 H > 5 T, but not in R3 and R4. Further analysis reveals that these single-frequency 
oscillations came from a relatively large Fermi extremal surface (Sf 2 = 49 T and 36 T for 
R1 and R2, respectively, see more details in SI). Given that these oscillations only appear 
in the samples with large carrier densities, we expect that they most likely come from a 
3 D electron band with a conventional energy dispersion. Further measurement in higher 
magnetic fields is prepared. 

At first glance, the longitudinal magnetoresistance (LMR) of R2 and R4 for H //a shows 
different profiles at low temperatures: R2 shows strong SdH oscillations up to 9 T, while 
R4 shows weak oscillations when H < 3.2 T. Beyond this field corresponding to the lowest 
quantized LL, R4’s LMR dramatically decreases with increasing H, ie ., R4 shows negative 
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LMR above the QL. The Lifshitz-Onsager rule is employed to understand the quantization 
of the Fermi surface cross-section area Sf as a function of the magnetic field B: 

= 27r ( n + 7) (!) 

where h is the Planck’s constant; e is the elementary charge; n is the LL index, and 7 
is the Onsager phase. The SdH oscillations of R2 originate from a small Fermi surface 
with Sp\ = 8.1 T, and the peaks at ± 6.8 T correspond to n — 1 (See Figj3] D and E). 
It is noteworthy that when T = 40 K, the oscillations fade out but the negative LMR still 
persists. On the other hand, the SdH oscillations in R4 originate from an even smaller Fermi 
surface with Sfi = 3.8 T, and a held larger than 3.2 T drives the system beyond the QL 
(See Figj3] F and G). The oscillations fade out above 10 K, but the negative LMR persists 
until 25 K, and then it is dramatically suppressed to a very small, positive signal above 
50 K. Although R1 and R3 have different Sf, their LMR manifests a similar behavior in 
rescaled H with respect to Sf- i), the negative LMR is always found to occur as long as 
the magnetic held is higher than the QL; ii), the negative LMR is not a part of oscillations 
since it survives at higher temperatures which destroy the oscillations. These phenomena 
described above have been observed in a dozen of samples with various sizes and n (see SI 
and Table HI for the samples R1 and R3). Based on the analysis we believe that the negative 
LMR is due to a band effect above the QL. It is noteworthy that the LMR of R3 and R4 
shows a shoulder-like anomaly around 6 T (indicated by the arrows in Fig. [2] and SI). These 
anomalies are located at the fields that are very close to the n = 1/2 in the LL indices 
(Fig. [4]C and E) for various samples with different Sf- It is likely due to the magnetic held 
induced spontaneous symmetry breaking which drives the system into a fractional filled LL. 
Further studies are needed to clarify this anomaly. 

In order to address the Fermi surface topology and the nature of the electrons contributing 
to the SdH oscillations, we measured the MR while rotating the direction of the magnetic 
held with respect to the three principal axes at 2 K (Fig. [3] A, B, C and SI). When the 
direction of the held changes from a to b or to c, the MR of R2 shows a very similar prohle 
with respect to the change of 0 and 0: the oscillations persist until 0 and 0 = 60°, while 
the negative MR is fully compensated by the linear contribution from the transversal MR 
for 0 and 0 > 30°. Figure [3] D and E show the resistivity ( pn ) versus the reciprocal of the 
held component along the a direction for R2. The oscillations’ period only depends on the 
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field component along the a direction, whenever the held is tilted to b or c. A similar angle 
dependent behaviors for R4 is shown in Fig. [3]F, G and SI. 

Based on the analysis of the SdH oscillations along different directions in R2, the fre¬ 
quencies are plotted as a function of the angles in Fig. [I] A and B. A clear 1 /cos(f> and 
l/cosip dependence up to 60° is observed on the Fermi surface Spi, which indicates that the 
electrons are distributed highly two-dimensionally and S F i has the geometry of a cylinder 
or a highly anisotropic ellipsoid. The LL indices for all four samples are shown in FigQ]C. 


The resistivity associated with the SdH oscillations is expressed as (ll2Qj: 


p H = p 0 [l + A(B,T)cos2tt(S f /B + -f)}, 


where 


A(B,T ) oc exp(—2ir 2 k B T D /ftujc) 


2n 2 k B T/hw c 


( 2 ) 


( 3 ) 


sinh(2TT 2 k B T / hwc) 

In Eq. [3l T B is the Dingle temperature; k B is the Boltzmann’s constant; and the cyclotron 
frequency u>c = eB/m*. Based on Eq. [2j the peak and valley positions of p zz are indicated 
as integer and half-integer, respectively. The interpolations of the lines n versus l/poH for 
all four samples show that the intercepts are clustered around 7 = —0.15 ± 0.05, although 
S F \ varies from 13.3 T to 3.8 T. This near-to-zero phase shift indicates a non-trivial Berry 
phase associated with Dirac or Weyl fermions ( 21 ], Q], The effective mass (m* = eB/ujc ) 
was obtained from the temperature dependence of the peak amplitude at n — 1 for R2 
and R4, yielding 0.075 and 0.064m e , respectively (Figj4]D). Other parameters including the 
Fermi wave vector (k F = yj -), the Fermi velocity (v F = fikp/m *) and the Fermi energy 
(Ep = v 2 F m *) are listed in Table HI These results, together with the Hall effect measurements, 
are consistent with the picture in which the Fermi levels of R3 and R4 lay at lower positions 
than those of R1 and R2 in a linearly dispersive electron band. 

According to our measurements on various samples, we believe that the negative LMR is 
not a current jetting effect due to highly anisotropic or inhomogeneous conductance 3, [231 


or improper contacts’ positions (see more details in SI), nor a weak anti-localization effect 
of impurity scattering which is usually observed in 2DEG in low magnetic fields [24]. Under 
the condition of piB < 1, the classical scattering effect generally changes little with respect 
to the magnetic held j^j]. The strong temperature dependence and the correlation with the 
QL indicate that such > 20% negative LMR must be cauesed by the band effect of Dirac 
fermions. Recently a weaker negative LMR was observed in Bi-Sb alloy with no indication 
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of SdH oscillations 


25]. Besides that, negative LMR has never been reported in other Dirac 


electron systems such as Bi and Cd : jAs 2 


3 - 29 ], 


The crystal structure of /3- Ag 2 Se differs from the other Dirac semimetals due to the 
absence of the central symmetry, which can lead to symmetry protected Weyl nodes induced 
by SOC Q- I n order to understand the chirality of electrons in Ag 2 Se, we analyzed the low- 
energy effective Hamiltonian from the space group P2i2i2i. In the absence of SOC, the single 
group of P2i2i2i for the spinless system possesses one-dimensional representation along T-X 
(X) and T-Y (A) axes in the Brillouin zone (BZ) 31]. The band structure possesses two 
Dirac-like bands around on the high-symmetric X and A lines, respectively (Fig. EDA). 
This result is consistent with the primary band structure calculation based on the density 
functional formalism (See SI). The bands hybridize each other and open gaps on the X and 
A lines in the presence of SOC (Fig. [5] B). One Dirac point at a high-symmetry line (X 
or A) splits to two Weyl points at different quadrants of the k x — k y plane. The 3D Weyl 
Hamiltonian near Ep characterizes the band dispersion and the spin texture in vicinity of 
the gapless point k 0 = (k x o, k y0 , 0) (Fig. |5]B). The other Weyl point near the A line in this 
quadrant has a higher energy and stays far from Ep. The two half-way screw operations of 
the Weyl point at k 0 yield two opposite chiral Weyl points at (—k x0 , k y 0 , 0) and (k x0 , —k y 0 , 0), 
while a time reversal operation gives the Weyl point with identical chirality at —ko (see Fig. 
0C). 

The strong SdH oscillations with a large anisotropy in Ag 2 Se indicate a large LL spacing 
when H //a (See SI). Above the QL, the electrons in a Weyl system entirely occupy the zeroth 
LL and only possess a momentum dispersion parallel to H 32,|33|. Such a chiral anomaly 
pumps electrons between Weyl points with the opposite chiralities at a rate proportional 
to the scalar product of the electric and the magnetic held E ■ B , giving rise to a negative 
LMR along the direction that connects the pair of Weyl points 34], [35]. According to our 
calculations, the LLs from opposite Weyl points in Ag 2 Se coexist around the projection 
point ±k 0 when H// a (Fig. ED D, see SI for details). Above the QL, the magnetic held 
further gives rise to a 26 shift along the direction of the held due to a dynamical chiral 
shift induced by the electron-electron interactions (Fig. [5] D) 36]. The large momentum 
mismatch between the Fermi points in the zeroth LL impedes the backscattering process. 
Because the density of states is proportional to the magnetic held in the zeroth LL in an one 
dimensional momentum space, the longitudinal magneto-conductance (cr = 1 /p zz ) increases 
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TABLE I: Parameters in samples Rl to R4. Spi, kpi and Sp 2 , kp 2 are determined from the 
SdH periods on a direction and b directions, respectively. Ep and vp are obtained from the T 
dependent amplitude of oscillations on a direction. Dash entries mean quantities not measured. 
N/A means quantities not observed. 



n 

cm -3 

S F 1 

T 

kpi 

A- 1 

7 

m* 

m e 

Ep 

meV 

Vp 

10 5 ms -1 

Sf2 

T 

kp2 

A- 1 

Rl 

4.9 x 10 19 

13.3 

0.02 

-0.14 ±0.04 

0.1 

30 

2.3 

49 

0.039 

R2 

- 

8.1 

0.016 

-0.18 ±0.07 

0.075 

27 

2.5 

36 

0.033 

R3 

- 

4.7 

0.012 

-0.17 ±0.01 

- 

- 

- 

N/A 

N/A 

R4 

2.4 x 10 18 

3.8 

0.011 

-0.18 ±0.01 

0.064 

18 

2.2 

N/A 

N/A 


linearly with /x 0 H above the QL (Fig. 0]E). This linear positive magneto-conductance in 
Ag 2 Se is not like an H 2 dependence due to the impurity scattering effect 
a characteristic of anomalous transport properties for a Weyl semimetal 


37 

32 


but more like 


34J. 


In summary, we synthesized macro-size single crystals of /3-phase Ag 2 Se and measured 
their resistance in magnetic fields. While the linear transversal MR was confirmed, the 
LMR showed strong SdH oscillations and a large, negative component beyond the QL. This 
anomalous MR was understood as an electron pumping effect induced by chiral anomaly in 
the zeroth LL. Our measurements also indicated that the negative LMR will be saturated at 
the base temperature when the magnetic field is far beyond the QL. Some first-step measure¬ 
ments hinted that the sample will become an insulator in a sufficiently high magnetic field. 
Finally, we infer that the vapour transfer growth method bears more efforts for synthesizing 
the single crystals of new materials, especially for the thermally unstable compounds. 


Methods 

Polycrystalline /3-Ag 2 Se was prepared from stoichiometric Ag and Se in a sealed quartz 
tube at 800 °C for two days. The yield was then ground and sealed in a long quartz 
ampoule in vacuum, and placed in a tube furnace with temperature gradient from 500 °C 
to a temperature close to the room temperature for 7 days. The single crystals yielded are 
ribbon-like and their sizes range from 20 /im x 20 [im x 2 mm to 0.5 mm x 0.5mm x 10mm 















(Fig. H] inset). The details of the crystal growth are discussed in SI part. Single-crystal X- 
ray diffraction data was collected on a Agilent SuperNova diffractometer with a micro-focus 
Mo source (A = 0.7107A). The data reduction was precessed by the CrysAlis Pro package. 
Multi-scan absorption correction was applied. All physical property characterizations were 


performed on a Quantum Design Physical Property Measurement System (PPMS-9) . 
first-principles calculations are based on the generalized gradient approximation (GGA) 


’he 


using the full-potential projected augmented wave method 
VASP package 


41 


39, 


40] as implemented in the 


43]. Electronic structure calculations were performed over 12 x 6 x 6 


Monkhorst-Pack 7-mesh included self-consistently. 
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FIG. 1: Temperature dependent resistivity of single crystalline f3- Ag2Se. The resistivity 
along a direction is similar as previously reported polycrystalline samples with similar carrier 
density n. Upper inset: the unit cell and the definition of angles 6, (j) and i/> among three principal 
axes. Central inset: an image of a crystal under microscope. The scale is 1 x 1 mm. 
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FIG. 2: MR (Ap/po) for the samples R2 and R4 along three principal axes at represen¬ 
tative temperatures. A p xx , A p yy and A p zz are for H// b, c and a, respectively. Since the raw 
data are highly symmetric, the treatment of pn = {p+H + P-h)/ 2 was not applied. The arrows in 
the lower right panel show the shoulder-like anomaly above the QL. 


14 
















0.0 0.2 0.4 0.6 0.8 1.0 

l/H 0 Hcos(j> (T 1 ) 


FIG. 3: Anisotropy of MR for the samples R2 and R4 at 2 K. Panel A, B and C: the MR 

for R2 when the direction of the field change among three principal axes at 2 K. No symmetry 
treatment of ±H was applied. The MR slightly changes when 6 changes from 0° to 90°. Panel D 
and E: the resistivity versus l/poH a = 1/ {/jlqH coscj)) and l/(p-oHcosip) for R2 at 2 K, respectively. 
Panel F and G: pu and —cPpu/dH 2 versus l/poH a = l/(poHcoscp) for R4 at 2 K. The red dashed 
line on the left represents the position of the shoulder-like anomaly beyond the QL (see more details 
in the main text). 
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FIG. 4: Analysis of SdH oscillations and negative LMR in Ag 2 Se. Panel A and B: S^i for 

R2 varying with cp and ip as l/coscp and l/cosip (red curve), respectively. Panel C: Landau level 
indice for all 4 samples. Two open symbols present the anomaly at n = 1/2 for R3 and R4. Panel 
D: Temperature dependence of the SdH oscillations amplitude at n = 1 for R2 and R4. Panel 
E: magneto-conductance for three samples showing negative LMR below 9 T at 2 K. The arrows 
indicate that a linear magneto-conductance occurs when the field is higher than the position of the 
shoulder-like anomaly. 
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FIG. 5: Schematic plots of energy bands of Ag 2 Se. Panel A and B: The energy bands for 
spinless and spinful system, respectively. Panel C: Weyl nodes in the BZ. Panel D: Landau levels 
for B // X when the dynamical chiral shift 2 6 is counted. 


17 















